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feedback controller can steer a pro- 
cess variable toward the desired set- 
point only if it can somehow pr edict 
the future effects of its current con- 
ktrol efforts. A model-based controller 
does so with the help of a mathematical repre- 
sentation of the process's behavior. A well-tuned 
PID loop uses an implicit model characterized 
by the values of the con- 
troller s P (proportional), 
I imtegrall, and D (deriva- 
tive! parameters. 

These two techniques 
compute their control 
efforts differently, but they 
both rely on the linearity of 
the process to anticipate 
how it is going to respond. 
A process is said to be lin- 
ear if the process variable 
increases by a factor of u 
when the control effort 
is increased by the same 
amount. And if two sepa- 
rate sequences of control 
efforts are added together 
and applied to a linear pro- 
cess, the resulting values 
of the process variable will 
always equal the sum of 
the values that would have 
resulted had the two con- 
trol efforts been applied 
separately. 



seconds has applied a different sequence of con- 
trol efforts to the same linear process. 

In case A, the controller has applied a single 
impulse with a magnitude of 1 unit (percent, 
degree, PS1, whatever) and a width of one cycle 
time (At seconds). The resulting fluctuations in 
the process variable are known as the process's 
impulse response or. in this particular case, its 



Superposition Principle 



Process response 



» A sluggish process will demonstrate a relatively slow unity im pulse response when the controller applies 
9 n on/off control effort with a magnitude of 1 (unity j over a period ot one cycle time (At seconds! . 



Control effort 



I Since the process is linear, it demonstrates the same impulse response magnified by a factorof u when 
the ointjoller applies a second impulse 11 times as large as the In st lnrhise«ample,u=3. 



Process response 



F 



Superposition 

This predictability gives 
rise to the Superposition 
Principle which governs 
the behavior of all linear 
processes. The "Superpo- 
sition Principle' graphic 
shows how it works in 
four situations where a 
computer-based control- 
ler with a cycle time of At 



At 

C When two impulses are applied to the process, the process's net response equals Hie sum of the two 
individual im pulse responses added together point by point. 



»<°) u{1). 



Hl- 
At 

Similarly, a sequence ot back-to-teck impulse 
At seoruids yields 3 sequence of super impos* 
achieved it the controller simply switches its 1 



Process response 



/ 



! *vith magnitudes u(0), uit). tt<2), ...applied every 
d impulse responses. Note that the same effect can be 
C'titrol effort directly to anew value at the end ot each 



eye to. ft need not turn its output all the way off then bscd on again between adjacent Impulses 
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Predict the Future 



1,0,0. ... u, 0,0, ... 



KO),Mt),u(23,. 



unity ftnpwtr response. 

The process in this exam- 
ple happens to be somewhat 
sluggish, so its unity impulse 
response rises and falls rela- 
tively slowly as the effects 
of the impulse wear off. This 
could represent any numlx-r 
of industrial processes, such as 
the temperature in a vat after 
a heating element has been 
turned on then off again, or the 
flow rate in a pipe after a valve 
has been opened then closed. 

Case B shows how increas- 
ing the magnitude of the 
impulse increases the magni- 
tude of the Impulse response 
but not its general shape. The 
second impulse is three times 
as large as the first, so the mag- 
nitude of the impulse response 
has been tripled. 

In ease C, both impulses 
have been applied to the pro- 
cess but at different times. Hie 
process's net response after the second impulse 
equals the sum of the two impulse responses 
added together point by point. The second 
impulse response has been effectively "superim- 
posed* on the first, hence the name of the prin- 
ciple that describes this phenomenon. 

Case D shows that a contiguous sequence of 
impulses with magnitudes of ujo), ujl), u(2l 
... applied to the process at times 0, At, 2At, ... 
has the same additive effect Each new impulse 
response adds to the impulse responses already 
in progress, and the magnitude of each is deter- 
mined by the magnitude of the impulse that 
caused it. The process's net response at any time 
is the sum of all the impulse responses that have 
been initiated up to that point 

Equivalent calculation 

Thanks to the Superposition Principle, a control- 
ler can predict how a linear process will respond 
to any sequence of control efforts, not just 
impulses. It also gives an algorithm for comput- 
ing the resulting values of the process variable 
as shown in Calculating the Process Response.' 



Calculating the Process Response 



A Here, a linear process has been eicited by a unity impulse ia single ccnlicl effort of 1 iinitj at time 0. 
The resulting unity impulse response Iras been recorded a: a sequence ol process variable 
measurements hrOj. h(1j, h(2., .. taken at times 0, it, 241. ... 



U,0,0... 



u-h(0),u-h(i).u-h(2).- 



I A second impulse u times as strong as trie first has magnified each value of Hie impulse response 
by a factor of u. 



h(0),...h(»), I" rK10)]J t*i«|,_. 



C Here the controller has applied two impulses - the first with a magnitude of 1 at time 0, the setond with 
a magnitude of u at time 1 0at The resulting process response equals the second Im pu Ise response 
superimposed on the first starling at time 1041. 



• J 



♦uU,-h/o> 



« In this more general case, the controller has applied a sequence of control efforts n(OI. 0(1), u(2), ... 
equivalent to an Impulse with a magnitude of ujO) at time 0, another vdlh a magnitude of u( 1 ) at time At. 
a third «Mh a magnitude of ui2) at time 241 etc. Each Impulse Initiates another Impulse response with 
a magnitude proportiona l to the magnitude of the corresponding contr ol effort. The net process 
response at any time equals the sum ol all the impulse responses initiated up to that point. 

fitmroe C<*ttrol Engineering 



This graphic depicts the same four situations, 
except that the control efforts and the corre- 
sponding process responses are represented by 
their numerical values rather than trend charts. 
Each data stream has been sampled and record- 
ed once every At seconds, hence the expression 
sampling interval often used to describe the con- 
troller's cycle time it. 

Case D shows the calculations required to 
compute the values of the process variable y(0), 
yit|, y(2), ... that would result from an arhitrary 
sequence of control efforts u(0), ujl), u(2), ... 
Specifically 



y|0) = ujOMi(0) 
y(lJ-u(Ol-h(ll-t 
vj2) = u.0l-ht2)H 



uill-hiOl 



etc. Each calculation gets successively longer as 
more and more impulses figure into the result. 
Fortunately, there's a convenient way to organize 
all these multiplication and addition operations, 
as shown in the 'Convolution* table, where two 
infinitely long 'numbers' 
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Convolution 





Column 1 


Column 2 


Column 3 


Column 4 




h(0) 


m 


h(2) 


h<3) 




u(0) 


u(1) 


.1(2) 


u(3) 


Row 1 


u<0) h(0) 


u(0)h<1) 


u(0) h(2) 


u(0) h(3) 


Row 2 




u(1)h(0) 


u(1>h(1) 


u(1)h(2) 


Row 3 






u(2)h(0) 


u(2)h(1) 


Row 4 








u(3)h(0) 



y(3) 



This table organizes the calculations required to compute the values of the process variable 
y(0). y(1j, y(2). ...that will result when control efforts mOi.Xti. i(2j. ... aie applied to a linear 
process Ml a unity Impulse response of h(0) . kj1 j. h(2). ... Row I is calculated by multiplying 
each value of the unity impulse response try nO). This represents the impulse response initi- 
ated by the first control effort wCn. Row 2 is calculated the same way except that the results 
are recorded starting in column 2 rather than column 1. and the multiplier is ui 1 ) This rep- 
resents the second impulse response initiated by control effort u!1 \ Each subsequent row is 
calculated similarly, starting one column further to the nght and using the next control effort as 
the multiplier. Adding the entries in each column from row 1 down io row k+1 sums all of the 
impulse responses Initiated from time through time kit and thereby yields the values of yiui 
through yik). The value ofyiki can also be calculated individually using 

y(k) = u(0) h(k)+u(1) h<k-1>+ +u(k-1) hOHulk) h<0) 
For example, the table's fourth column show that 

y(3)= u(0) h(3Hu(1 ) h(2Hu(2) tK1)+u(3) h(O) 

Source: Ccotrv' Engineering 



B = NO), l>t2), ... 
and 

V = u(O), U(ll, u(2), ... 

are multiplied' together to compute 

V = yjO), yil), y(2), ... 

using the familiar long multiplication algorithm 
but with data points h|OI. h(ll, h(2!, ... and 11(01, 
iilll, ti(2l. .. instead of individual digits- This 
calculation, known as conwlution. is actually the 
mirror image of long multiplication. The multi- 
plication and addition steps are the same, but it 
does not involve any carry-over from one column 
to the next. It is typically written as Y = JI*U 
where "*" is the convolution operator 

Convolution is the basis for an entire math- 
ematical discipline known as linear systems analy- 
sis. It gives control engineers a powerful tool for 
analyzing the behavior of linear processes and 
designing feedback controllers that can predict 
the future. 

Vance Van Doren, Ph.D., P.E., is senior 
editor for Control Engineering. He can be reached 
at vanceiS-control.com. 
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Analyzing 

Closed-Loop Behavior 
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The second in a two-part series explaining how linear process behavior 
can be predicted mathematically. 
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The repeatable behavior of a linear 
process allows a process controller to 
predict the future effects of its cur- 
rent efforts. As described in 'Process 
Controllers Predict the Future |March 
20081, a controller's predictive abilities can be 
reduced to a mathematical formula by conduct- 
ing the following experiment: 

1) Stimulate the process in open-loop mode 
by forcing the controller to apply a unity impuke 
that is a single control effort that is one unit in 
magnitude and At seconds in duration. 

2) Determine the process s open-loop impulse 
response by measuring the process variable every 
At seconds after the impulse 

Analyzing feedback loops 

S = s(0). s(1). tO). ... A E » <H0>, <K1 ). e<2). ... _ 



31 Label those measurements in chronologi- 
cal order as h(0) h(1) h(2) etc and construct 
a single infinitely long 'number' H using those 
values as "digits' That is, let 

H = h(0). h(1). h(2). ... 

4| Start over and let the controller apply an 
arbitrary sequence of control efforts to the pro- 
cess. Label those values u(0) u(1) u(2) etc.. and 
use them to construct a second infinitely long 
number U: 

U = u(0). u(1). u(2), ... 




Controller 

a 



U = U(0).U(1).U<2).. 



Setpoint 
Error 
Control effort 

Disturbances 

Process variable 



Y = yi0).y(1),yt2),... 



D»d(0),d(l), d<2}, ... 



The behavior ot this closed-loop control system can be predicted by constructing tour -numbers' -S, 0. G. andH - with 
TMgits' consisting of sampled values ot the setpoint. the disturbance, the open-loop impulse response ot the process, and 
the open-loop impulse response of the conlroller. respectively. Thesequenceol process variable measurements Y that will 
result from any given combination o/S. D. G, and H can then be computed from 



Frequency Domain Analysis 
Explained." October 2005 



G-H * s + 



1+G*H 1+G'H 
vihere mestar c.i indicates convolution and the two quotients indicate deconwlution. 
Sou/ca. Ccntrol Engineering 
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Stability test 



■ mini. 



The histogram shews each of the digits calculated tor the transfer function 
T G'H 



1 + G'H 



that corresponds to one particular closed-loop system composed of a controler with an 
open-loop impulse response of G and a process with an open-loop impulse response of H. 
The graph of T, shows how the process variable would react if the closed-loop system 
were to be simulated with a single unity impulse (that is. it the setpoint sequence were to 
be changed to S ■ 1 , o, 0. ...). The shape of the graph wil vary for dtbrent combinations 
of processes and controllers, butiflhe closed-loop system Is stable, allot the digits of T t 
wil have trite values. 



1 1 n 1 1 1 « . 



■ 1 1 1 1 

iiiiinii. 




Hill 



I 



For an unstable closedloop system, the histogram of T, wil not turn out to be entirely fhte. 
and the system's cbsed-loop impulse response wil grow larger and larger over tme. 

Source Contrd Engineering 



5) Compute Y = H*U using convolution, which 
multiplies H and U together just as if they were 
normal multi-digit numbers. That is, use normal 
multiplication and addition to compute 

y(0) = u(0) ■ h(0) 

y(1) = u(0)h(1) + u(1)h(0) 

y(2) = u(0) • h(2) + u(1) • h(1) + u(2) • h(0) 

and so on until the kth sampling interval when 

y(k) = u(0) • h{k) + u(1) • h{k-1) + ... 
+u(k-1) ■ h<1) + u(k) • h(0) 

If the process is truly linear and there are no 
other influences affecting the process variable 
then these computed values should match the 



actual process variable measurements that will 
eventually result from the control efforts U. 

Conversely the equations in step 5 could 
be solved for the open-loop control efforts U 
that would be necessary to generate a speci- 
fied sequence of process variable measure- 
ments Y That is, pick the desired values of 

Y = y(0), y(l).y(2},... 

and let 



u ( o) = yi£) 
m 



u(1) = 



. 1 

h(0) 



[y(l)-u(0)h(i)3 



u(2) =J_ [y(2) -u(0) - h(2) - U(1) • h(1)J 
h(0) 

and so on until the kth sampling interval when 

u( k) = _L ty(k) - u(0) ■ h(k) - u<i) ■ h(k -1) - . . . 

m 

-u(K-l)-h(l)] 

This operation, known appropriately 
enough as deconvohition, is essentially the 
long division algorithm applied to dividing Y 
by H to get U = Y/H. See the "Deconvolution 
example" graphic. 

Closed-loop analysis 

In practice, process controllers typically com- 
pute their control efforts using algorithms 
more sophisticated than open-loop control, 
but convolution and deconvolution can be 
used to analyze the behavior of closed-loop 
control systems as well. 
Consider, for example, the basic feedback 
control system shown in the "Analyzing feed- 
back loops' graphic. The process variable is 
sampled at regular intervals and each measure- 
ment is subtracted from the setpoint to gener- 
ate a sequence of error measurements Those 
are fed into a controller that in turn generates 
a sequence of control efforts chosen to drive the 
process variable toward the setpoint 

The resulting process variable measurements 
Y will equal the results of the controller s efforts 
added to any uncontrollable disturbances D that 
might also be affecting the process. In terms of the 
convolution operator, that sum can be expressed as 

Y = H ' U + D 

If the controller is also a linear process land 
most are I, then 
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Deconvolution example 



denominator (B) 



16. 4. 


1, 




|1, 0, 


0, 


0, ... 


16. -4. 


0. 


0. ... 


4. 


0. 


0. ... 


4. 


-1. 


0. ... 





1. 


0, 




1, 


... 







ft 

i, ... 



<- quotient (A/B) 
<- numerator (A) 



U = G 1 E 

where G is the impulse response of the controller, S is 
the setpoint sequence and 

E = S-Y 

is the error sequence. These three relationships can be 
combined by means of simple algebra treating convo- 
lution just like multiplication and deconvolution just 
like division. That is, 

Y = G * H * (S - Y) + D 



This table shows how two sequences of sampled data values ,A = 16. 
0, 0. ... and B = 1, 0, 0, can be deconmlved to compute a third 
sequence ( A/B = 16, 4, 1, ! i, Each row is generated using tlie long divi- 
sion algorithn except that individual "digits' can be fractions greater than 9 
and less than zero, and there is no carry-over from column to column. R>r 
example, the fifth row is 4 times the denominator, and the sixth row is the 
difference between the fourth and fifth rows The midtipher of 4 is chosen to 
produce a zero in the first column of the sixth row, and tliat 4 appears in the 
se.:ond column of the quotient. 



G* H 
1 + G * H 



S+_ 



1 +G*H 
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This formula predicts the process variable measure- 
ments Y that will result from any given combination of 
setpo'mts S and disturbances D when a process with 
an open-loop impulse response of H is controlled 
by a feedback controller with an open-loop impulse 
response of G. It is the fundamental mathematical tool 
used by control theorists to design controllers and ana- 
lyze their closed-loop behavior. 
For example, the quantity 



1 + G * H 



known as the transfer function between the setpoint and the pro- 
cess variable can be used to determine how the process variable 
Y will react to any given sequence of setpoints S in the absence 
of disturbances. That is, let D=0 so that 



The sequence T s is also the impulse response of the closed- 
loop system. That is, T, shows what the process variable would 
look like if a unity impulse were to be substituted for the set- 
point while the controller is on-line controlling the process 

This phenomenon yields a test for closed-loop stability of a 
feedback control system. If the impulse responses of the con- 
troller and open-loop process are both known then their con- 
volution product G*H can be used to compute T,. If all of the 
digits of T, turn out to be finite then the closed-loop system 
will be stable. Otherwise, Y - T S "S will grow forever no matter 
what values are chosen for S and the closed-loop system will 
be unstable. See the 'Stability test" graphic. 

On the other hand, if the setpoint S is zero and the controller 
is only concerned with counteracting the effects of disturbances 
then 

Y = T ' D 

where the quantity 
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1 + G " H 

is the transfer function between the dis- 
turbance sequence D and the process vari- 
able Y It shows what the process variable 
would look like if the closed-loop system 
were to be stimulated with a single unity- 
impulse disturbance It too can be comput- 



ed from G'H though it will typically have 
a different value than T . This is why a 
closed-loop system s disturbance response 
typically differs from its setpoint response 
though if one is stable so will the other. 

The transfer functions T, and T D can 
also be used to design the controller if H is 
known and G has not already been select- 



ed For example if the controller will be 
used primarily to track setpoint changes 
T s can be assigned a suitable value that 
will produce some desirable closed-loop 
setpoint response The controller required 
to achieve that particular setpoint response 
would then be given by 
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H*(1"T S ) 

Similarly, if the controller will be used 
primarily to reject disturbances, then T D 
can be assigned a suitable value, and the 
contr oller s impulse response can be set to 



H'T„ 



Devising a controller that will dem- 
onstrate either of these impulse response 
sequences exactly can be I bit tricky but 
there are mathematical techniques collec- 
tively known as transform theory that can 
help. The first step is to express G and H as 
fractions with numerators and denomina- 
tors of finite length For example, the infi- 
nitely long open-loop impulse response 

H = 16, 4. 1, %, ... 
can also be expressed as H=A/B with 
A= 1S 
B = 1,-% 

in much the same way as the infinitely long 
decimal value 0.14285... can be expressed 
as 1/7, In fact, deconvolution and normal 
division are essentially the same algorithm 
See the "Deconvolution example' graphic. 

The simplified version of H can also 
be converted into a Z transform that can 
be used to analyze a closed-loop system s 
behavior without having to perform any 
convolution operations involving infi- 
nitely long sequences of sampled data 
Unfortunately it takes considerably more 
mathematical expertise to understand 
how and why Z transforms work. 

Z transforms can be converted into 
Laplace and Founer transforms by decreas- 
ing the sampling interval At to zero These 
transforms can be used to analyze the 
behavior of analog control systems where 
data is processed in a continuous stream 
rather than at regular intervals e» 

Vance VanDoren lcontroleng0ansn.com! 
is consulting editor to Control Engineering. 
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